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ABSTRACT 

Using the Mid-Infrared Spectrometer (MIRS) on board the Infrared Telescope in 
Space (IRTS) we obtained the 4.5 to 11.7 /jm spectra of the stellar populations and 
diffuse interstellar medium in the Galactic bulge (/ f» 8.7°, b « 2.9, 4.0, 4.7, and 5.7°). 
Below galactic latitudes of 4.0°, the mid-infrared background spectra in the bulge are 
similar to the spectra of M and K giants. The UIR emission bands (6.2, 7.7, 8.6, and 
11.3 /um) are also detected in these regions and likely arise from the diffuse interstellar 
medium in the disk. Above galactic latitudes of 4.0°, the mid-infrared background 
spectra are similar to the spectra of those oxygen-rich evolved stars with high mass-loss 
rates detected by IRAS. One likely interpretation is that this background emission 
arises predominantly from those stars with very low luminosities that have not been 
detected by IRAS. The age for such low-luminosity evolved stars could be 15 Gyr, and 
the existence of a large number of evolved stars with high mass-loss rates in the bulge 
has a significant impact on our understanding of the stellar content in the Galactic 
bulge. 
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1. Introduction 



The Galactic bulge is characterized by a diversity of stellar populations: main-sequence stars, 
K & M giants, Asymptotic Giant Branch (AGB) stars, planetary nebulae and white dwarfs. Bulge 
ages derived from these different stellar populations have big discrepancies. Studies of the RR 
Lyrae variables, the main-sequence turnoff, the luminosity of M giants, the mass of planetary 
nebulae, the characteristics of Mira variables, and AGB stars with large mass-loss rates give bulge 
ages ranging from 5 to 15 Gyr (e.g. Gratton et al. 1986, Terndrup 1988, Frogel & Whitford 1987, 
Kinman et al. 1988, Feast & Whitelock 1987, and van der Veen & Habing 1990, 1 Gyr = 10 9 
yr). Knowing the stellar populations and their properties in the Galactic bulge would not only 
provide useful information on the bulge formation, but would also help us to construct stellar 
synthesis models of spheroidal galaxies, for which only the integrated light can be observed. From 
the observations of the Diffuse Infrared Background Experiment (DIRBE) on board the Cosmic 
Background Explorer (COBE), Arendt et al. (1994) found that the near- infrared background 
radiation in the Galactic bulge is similar to that of late-K and M giants. Infrared Astronomical 
Satellite (IRAS) detected a large number of mid-infrared bright stars in the Galactic bulge 
(Habing et al. 1985), however, the IRAS did not have the sensitivity to detect those mid-infrared 
faint stars that could dominate the mid-infrared background radiation in the Galactic bulge. 
Furthermore, the IRAS Sky Survey Atlas (ISSA: Wheelock et al. 1994) flux densities at 12 and 25 
fim cannot provide unambiguous color data to study the stellar populations, which are likely to 
be contaminated by the interstellar unidentified infrared (UIR) bands and small grains emission, 
respectively. Thus, spectrophotometric or high spatial-resolution photometric observations in the 
mid-infrared band (e.g. 4-25 fim) are required for the study of the mid-infrared background 
radiation in the Galactic bulge. In this paper we present results from our mid-infrared (4.5 to 11.7 
lim) spectrophotometric observations in parts of the Galactic bulge, and we also discuss the nature 
of the stellar population there that dominates the mid-infrared background radiation. Throughout 
our paper we assume a Sun-to-Galactic center distance of 8.5 kpc. 



2. Observations and Results 

The Mid-Infrared Spectrometer (MIRS) was one of the four science focal-plane instruments 
that flew aboard the orbiting Infrared Telescope in Space (IRTS). This telescope was a joint 
Japanese Space Agency (ISAS)/NASA project that was launched on March 18, 1995 and surveyed 
~ 7% of the sky over the course of its 26 day mission life (Murakami et al. 1996). The MIRS had 
a beam size of 8' x 8' and operated over a wavelength range of 4.5 to 11.7 fim with a resolution 
of 0.23 to 0.36 /im (Roellig et al. 1994). During the course of its mission, the MIRS observed 
parts of the Galactic bulge and we present here the spectra at I « 8.7°, b rs 2.9, 4.0, 4.7, and 
5.7°. The absolute calibration is about 10% in this preliminary phase (see Tanabe et al. 1996 
for details of the MIRS calibration). Selected regions were chosen for this study which did not 
include bright IRAS point sources and the zodiacal background emission has been subtracted. 
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To subtract the zodiacal background, we first selected zodiacal background spectra that were 
observed by the MIRS at regions (far away from the galactic plane) that have the same ecliptic 
latitudes but different longitudes as those we chose in the bulge. We then corrected the brightness 
of these observed zodiacal background spectra to match the zodiacal brightness at the same solar 
elongations as those in the bulge by using a IRAS zodiacal cloud model (Wheelock et al. 1994). 
Finally, these model zodiacal background spectra subtracted from our selected bulge spectra 
(bulge + zodiacal emission). The zodiacal cloud emission at 12 fim in our selected regions is 
between 5 and 9 times higher than the bulge emission there. 

We estimated the extinctions to our observed regions for the purpose of checking whether 
extinction corrections on our spectra are necessary. Based on the DIRBE/COBE data, Arendt 
et al. (1994) found that the unreddened near-infrared colors in the Galactic bulge are similar to 
those of late-K and M giants. Thus, the optical depth at 1.25 fim can be estimated as 

_ ^[(/i,25 M m/^2.2^m)/0-95] 
Tl.25 M m " _ Q 

where I is the measured DIRBE intensity, 0.95 is the ratio of the flux densities at 1.25 and 
2.2 fim of the unreddened late-K and M giants, and -0.603 = A2. 2^^1.25^771 - 1 is the reddening 
between 1.25 and 2.2 fim and is adapted from the Rieke & Lebofsky (1985) interstellar extinction 
law. Based on equation (1) and the final DIRBE/COBE Mission-averaged Pass 3B data with 
zodiacal light removed as described by Kelsall et al. (1998), the T\.2^ m are calculated by assuming 
that the extinction lies in front of the emission sources along the line of sight. It should be note 
that this assumption minimizes the total line-of-sight extinction derived from the reddening. Table 
1 summaries the T\,2b^m as derived in the closest DIRBE pixels to our observed positions. The size 
of a DIRBE pixel is 0.32° x 0.32° and the size of DIRBE beam is 0.7° x 0.7°. With these derived 
T\.2bam values, the optical depths at 4.9 fj,m (r4.9„ m ) and 10 /jm silicate (T s ui ca t e ) are calculated 
from the Rieke & Lebofsky (1985) extinction law and are also summarized in Table 1. In Table 
1, the extinctions at /or around our observed positions are small (1 - 10% in the MIRS spectral 
range), therefore, extinction corrections on our spectra were not applied in the following analysis. 

Figure 1 shows the observed spectra of the Galactic bulge at our selected positions. In Figure 
1 (a) & (b), it can be seen that the mid-infrared background spectra in the bulge at Galactic 
latitudes of 2.9° and 4.0° are similar to those of M and K giants observed by MIRS, which do not 
have thick dust shells and the mid-infrared emission mostly arises from the stellar photospheres 
at temperatures of a few thousands K (Yamamura et al. 1996a, b). The UIR emission bands (6.2, 
7.7, 8.6, and 11.3 /im) are also detected in these regions and they likely arise from the diffuse 
interstellar medium in the disk (see e.g. Onaka et al. 1996 for the observations of interstellar 
UIR emission bands in the galactic plane). Figure 2 shows the spectrum of one of the M giants 
observed by MIRS and it can be seen that its spectral shape is similar to the stellar spectra in 
Figure 1 (a) &i (b). On the other hand, in Figure 1 (c) & (d) the mid-infrared background spectra 
at Galactic latitudes of 4.7° and 5.7° are similar to those of oxygen-rich AGB stars with thick 



(1) 
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dust shells, which have a typical silicate absorption feature at 10 (e.g. Forrest et al. 1978 and 
Herman et al. 1984). For comparison, Figure 3 shows the spectrum of an oxygen-rich AGB star 
with silicate absorption feature observed by MIRS. Even though the MIRS spectra in Figure 1 (c) 
& (d) are a bit noisy the silicate absorption feature is real. This is because the depression in the 8 
to 11 /jm range in both of the spectra is ~ -75% ±25%, which is larger than the 10% calibration 
uncertainty of the MIRS. Furthermore, these silicate absorption features should not arise from 
the line-of-sight extinction, because the optical depths of the interstellar silicate absorption at/or 
around our observed positions are much less than one (see Table 1). Also shown in Figure 1 are 
the IRAS Sky Survey Atlas (ISSA) 12 fim fluxes. In Figure 1 (a) & (b), the ISSA 12 fim flux 
seems to be mostly due to UIR emission bands at 7.7, 8.6, and 11.3 /im. On the other hand, the 
ISSA 12 /xm flux in Figure 1 (c) Sz (d) is mostly due to circumstellar dust emission. 



3. Discussion 

3.1. Evolved stars as seen by IRAS 

Habing (1987) (hereafter Habing87) has shown that the IRAS 25/12 /xm flux ratio for IRAS 
point sources in the Galactic bulge can be used to select those late-type stars with high mass-loss 
rates to produce an estimate of the bulge age. Habing87 assumed that the total integrated flux of 
these stars can be represented as 3zT\,(12 /mi), where v = 2.5 x 10 13 Hz and F u (12 fiva) is the 
measured IRAS 12 flux density. The factor 3 is defined as the infrared bolometric correction 
(see van der Veen Sz Breukers 1989). Habing87 proposed that these selected objects in the bulge 
region are Mira and OH/IR stars, and found that they have a luminosity function which peaks 
at ~ 4000 Lq. For an AGB star with a luminosity of 4000 L , Habing87 estimated that it has 
an initial mass of 1.7 M and age of ~ 1.3 Gyr (from zero-age main sequence to the first giant 
branch, Z = 0.04, and Y = 0.3) - much lower than that estimated from observations in Baade's 
windows (> 10 Gyr) (see e.g. Frogel 1988). Later study by van der Veen Sz Habing (1990), with 
more accurate selection criteria to select AGB stars from the IRAS data, find that the AGB stars 
with high mass-loss rates in the Galactic bulge have a luminosity function that peaks at ~ 5500 
Lq. However, with better knowledge of the evolution of AGB stars, they found that the initial 
mass of AGB stars with 5500 L is only 1.4 M and their age is 7 - 15 Gyr (Z = 0.04 and Y = 
0.2). Harmon Sz Gilmore (1987), however, have found that the AGB stars as seen by IRAS in the 
bulge have a luminosity function that peaks at ~ 6000 and 7600 L Q . They estimated these AGB 
stars have initial masses of 1.3 and 1.5 M and ages of 10 and 5 Gyr (Z = 0.1 and Y = 0.25), 
respectively. From the above IRAS studies, the ages of AGB stars in the Galactic bulge have a 
large range of 5 - 15 Gyr. The main reason for this large range is due to the high uncertainty of 
the mass-loss history of these AGB stars in the bulge. 
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3.2. Evolved stars as seen by IRTS 

In Figure 1 (c) & (d), the mid-infrared background spectra at Galactic latitudes of 4.7° and 
5.7° are similar to the spectra of evolved stars with circumstellar silicate absorption. Therefore, 
the stars in these regions that produce the mid-infrared background emission are very similar to 
those AGB stars with large mass- loss rates detected by IRAS. Since there are no bright IRAS 
12 /im point sources in these bulge regions that can account for the flux levels measured by the 
MIRS, one likely interpretation of this background emission is that it arises predominantly from 
AGB stars with large mass-loss rates, but with very low luminosities that could not have been 
detected by IRAS. From the IRAS Faint Source and Point Source Catalogs we find that there are 
less than 50 detected IRAS point sources per square degree at our observed positions in Figure 1 
(c) & (d), which is less than the threshold value (50 per square degree) of IRAS confusion limit at 
12 /jm (IRAS Explanatory Supplement: Beichman et al. 1988). Therefore, the IRAS observations 
at our observed positions in Figure 1 (c) & (d) are not confusion-limited. Away from confused 
regions of the sky the sensitivity of IRAS at 12 fj,m is 0.5 Jy (Beichman et al. 1988). For an 
upper-limit estimate of the stellar properties (e.g. luminosity and initial mass), we assume that 
each of these AGB stars has a 12 /j,m flux density of 0.5 Jy, and that this flux density corresponds 
to a luminosity of 850 Lq in the bulge. From the MIRS spectra in Figure 1 (c) & (d), the flux level 
of 30 Jy at 11.7 ^m implies that there are at least 60 AGB stars in the 8' x 8' MIRS beam (or ~ 
3375 stars deg -2 ), each averaging ~ 850 L Q (cf. Perault et al. 1996 of ISO survey observations 
of the Galactic plane at / « -45° and 6^0°, where they detected ~ 1500 M & K giants deg -2 
at 15 fim). Comparing the spectra in Figure 1, it can be seen that M & K giants dominate the 
mid-infrared emission in the lower latitudes of the bulge (Fig. la & b), while low luminosity AGB 
stars with large mass-loss rates dominate the mid-infrared emission in the higher latitudes (Fig. 
lc &: d). These differences indicate that the bulge has a gradient in the stellar populations such 
that the relative numbers of M & K giants to AGB stars decrease with increasing latitude. 

The relation between the maximum luminosity in the late-stage of AGB star and its initial 
mass were discussed by Iben & Renzini (1983), and can be written as 

L max ^ 59250 [0.53 t?- - 082 + 0.15 (M; - 1) ??" a35 - 0.495] , (2) 

where L max is the maximum luminosity (in unit of L ) reached along the AGB evolution, the 
parameter rj is proportional to the mass-loss rate of the star and has a value of unity with a factor 
of 3 uncertainty (1/3 < r/ < 3), and Mj is the initial mass (in units of Mq) of the AGB star. Here, 
we choose r\ = 1 for estimates of the initial mass of our detected AGB stars, and from equation (2) 
we find that these AGB stars have initial masses of 0.86 M . The age of these 0.86 M stars can 
be estimated if their metallicities are known, unfortunately however, our observations presented 
here do not provide any useful information about their metallicities. Just for a comparison with 
the past IRAS studies mentioned in the previous section, we assume that these AGB stars have 
metallicities of 1.9 times solar (Z = 0.03). From the stellar evolution model of VandenBerg & 
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Laskarides (1987), the age of a 0.86 M star with Z = 0.03 and Y = 0.35 is ~ 15 Gyr. 

In summary, the major results of this paper are the following: (1) there is a gradient in the 
stellar populations of the Galactic bulge, where the relative numbers of M & K giants to AGB 
stars decrease with increasing latitude; and (2) there could be a lot of (3375 per square degree) 
old age (~ 15 Gyr) and low luminosity (~ 850 L Q ) AGB stars with large mass-loss rates at the 
higher latitude field of the bulge. 

We are grateful to Rick Arendt for providing the optical depths at 1.25 summarized 
in Table 1 and many useful comments which helped in improving the paper. We thank Martin 
Cohen for providing the calibrated spectra of standard stars for the calibration of MIRS. We also 
thank IPAC for their help in the pointing reconstruction of the IRTS and the entire IRTS team for 
their efforts in ensuring the success of the IRTS mission. I. Y. is supported by the JSPS Research 
Fellowships for Young Scientists. The COBE data sets were developed by the NASA Goddard 
Space Flight Center under the guidance of the COBE Science Working Group and were provided 
by the NSSDC. 
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Table 1. The extinctions in the closest DIRBE pixels to our observed positions. 



galactic coordinates (/ , b) 


(8.78° , 2.71°) 


(8.70° , 4.16°) 


(8.37° , 4.72°) 


(8.62° , 5.60°) 


T 1.25^™ 


0.415 


0.247 


0.200 


0.122 


T~4.9/im 


0.034 


0.020 


0.016 


0.009 


Tsilicate 


0.096 


0.057 


0.046 


0.028 
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Fig. 1 — The observed MIRS spectra at I « 8.7°, b « 2.9, 4.0, 4.7, and 5.7°. The measured fluxes 
have been averaged in units of 5 second and this resulted in an effective beam size of 8' x 20'. Also 
shown in the figure, the * symbol, is the IRAS Sky Survey Atlas (ISSA) 12 /im flux. The bandpass 
of the IRAS 12 fj,m channel was 8 to 15 /im, and the resolution of the ISSA 12 /tm data is 4'. The 
total flux from the IRAS point sources contributes no more than 10% of the ISSA 12 /im flux in 
these regions. 

Fig. 2. — The spectrum of Z UMa (M5 giant) observed by MIRS and adapted from Yamamura et 
al. (1996b). 



Fig. 3.— The spectrum of OH 138.0 +7.2 (an OH/IR star) observed by MIRS. 
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